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ABSTRACT 

Images from the Hubble Space Telescope Advanced Camera for Surveys are used 
to carry out a new photometric study of the globular clusters (GCs) in M104, the 
Sombrero galaxy. The primary focus of our study is the characteristic distribution 
function of linear sizes (SDF) of the GCs. We measure the effective radii for 652 clusters 
with PSF-convolved King and Wilson dynamical model fits. The SDF is remarkably 
similar to those measured for other large galaxies of all types, adding strong support 
to the view that it is a "universal" feature of globular cluster systems. 

We use the Sombrero and Milky Way data and the formation models of Baumgardt 
& Kroupa (2007) to develop a more general interpretation of the size distribution 
function for globular clusters. We propose that the shape of the SDF that we see 
today for GCs is strongly influenced by the early rapid mass loss during their star 
forming stage, coupled with stochastic differences from cluster to cluster in the star 
formation efficiency (SFE) and their initial sizes. We find that the observed SDF shape 
can be accurately predicted by a simple model in which the protocluster clouds had 
characteristic sizes of 0.9 ± 0.1 pc and SFEs of 0.3 ± 0.07. 

The colors and luminosities of the M104 clusters show the clearly defined classic 
bimodal form. The blue sequence exhibits a mass/metallicity relation (MMR), follow- 
ing a scaling of heavy-element abundance with luminosity of Z ^ L'^-^ very similar to 
what has been found in most giant elliptical galaxies. A quantitative self-enrichment 
model provides a good first-order match to the data for the same initial SFE and 
protocluster size that were required to explain the SDF. 

We also discuss various forms of the globular cluster Fundamental Plane (FP) 
of structural parameters, and show that useful tests of it can be extended to galaxies 
beyond the Local Group. The M104 clusters strongly resemble those of the Milky Way 
and other nearby systems in terms of such test quantities as integrated surface density 
and binding energy. 

Key words: galaxies: star clusters - galaxies: individual (M104) - globular clusters: 
general 



1 INTRODUCTION 

Our knowledge of globular cluster (GC) systems in spi- 
ral and disk galaxies is far more limited than those in 
elliptical galaxies, and consists primarily of the samples 
in the Milky Way and in M3 1, along with a handful 



of more distant disk gala xies (iKissler-Pa tig et al. 19991: 



iGoudfrooii et al. |[2007 : Chandar et al. 112004 : .Rhode et al. I 



^ This work was based on observations with the NASA/ESA 
Hubble Space Telescope, obtained at the Space Telescope Science 
Institute, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under NASA contract NAS 5-26555. 



2007 : ISpitler et al. Il2006l:l Mora et al. 1l2007l : ICantiello et al. I 



20071 : iDeGraaff et al. I 120071 '). These should be compared 



with the many extensive studies of the GC systems in el- 
liptical galaxies, particularly the giant E's in which GCs 
can be found by the thousands (e.g. Larsen et al. ' '200l|; 
iBrodie fc Strader | [2006: Peng ct al. 2006; Harris 2009a, to 
name only a few). But if we are to probe the systematic 
properties of the GCs themselves and to understand fur- 
ther the common nature of these classically ancient objects, 
we need to investigate them closely across all types of host 
galaxies. 

The Sombrero galaxy (M104 = NGC 4594), a giant Sa 
system, is the disk galaxy with the largest known popula- 



2 Harris et al. 



tion of globular clusters (GCs) and for that reason alone 
is a particularly attractive system. The fact that it is dis- 
tinctly closer (d = 9 Mpc) than the Virgo cluster, and also 
nearly edge-on to our line of sight so that our view of its 
GC population is minimally affected by issues of reddening 
and contamination from the disk itself, m akes M104 an al- 
most unique target. iRhode fc Zepf] l|2004l ). from wide- field 
ground-based imaging, derive a total GC population of 1900 
clusters. They find that the spatial extent of the system is at 
least as large as 40 kpc i n projected galactoc entric distance. 

In an earlier study (jSpitler et al. Ibood which we refer 
to as Paper I), a first round of HST-based photometry of 
M104 GCs and measur e ments of their scale sizes was pre- 
sented. IChandar et al. I (|2007l ) used the same material to 
evaluate the trend of mean GC density with galactocentric 
distance. The raw data consisted of a special imaging dataset 
of six fields taken with the ACS Wide Field Channel as part 
of a Hubble Heritage Project on M104 (PL K. Noll, PID 
9714). This mosaic was observed in the standard BVR fil- 
ters (F435W, F555W, F625W) and covers a total field size 
of about 600" x 400" centered on M104. In Paper I, pho- 
tometry from the individual images was used to discuss the 
distribution of the GCs in luminosity and color, their spa- 
tial distributions around the galaxy, the differences between 
the red (metal-rich) and blue (metal-poor) GC subsystems, 
and the correlations of cluster sizes with GC luminosity and 
galactocentric distance. In addition. Paper I also provided 
the first indication that the metal-poor GCs in disk galaxies 
showed a correlation between mass and metallicity that had 
already been uncovered for elliptical galaxies. 

In the present study, we use the reprocessed mosaic of 
ACS images to measure the GCs in M104. A particular em- 
phasis of this study is the use of newer algorithms to take 
a deeper look into the intrinsic scale sizes of the GCs, and 
carry out a more extensive set of comparisons with other 
galaxies. This work leads us into a physical model for the 
origin of the GC size distribution that we observe today. 

In Section 2 we describe in detail the photometry and 
size measurements. In Section 3, we discuss the correlations 
of GC scale size with luminosity, metallicity, and galacto- 
centric distance, and show that versions of the "Fundamen- 
tal Plane" of structural parameters can be accurately con- 
structed for GCs in galaxies this distant. In Section 4, we 
introduce a quantitative model for the physical origin of the 
GC size distribution. The present-day GC size must be the 
result, not only of the slow dynamical evolution of the clus- 
ter in the tidal field of the galaxy, but also of the expan- 
sion of the protoclusters during their early phase of star 
formation and rapid gas expulsion. The key parameters in 
this model are the original protocluster size, the star for- 
mation efficiency, and (just as importantly) the dispersions 
in these two quantities. Finally, in Section 5 we redetermine 
the correlations of GC color (metallicity) with luminosity for 
both the bimodal sequences with our new photometry. We 
conclude more strongly than before that this galaxy has a 
mass/metallicity relation quite similar to those in the giant 
ellipticals. Section 6 gives a brief summary of our findings. 
(Readers interested primarily in the results may skip directly 
to Section 3 without loss of continuity.) 

Throughout this paper, we adopt the distance modu- 
lus (m - M)o = 29.77 (d = 9.0 Mpc; see Table 4 in Paper 
I) along with a foreground reddening Eb-v = 0.05 from 



the NASA/IPAC Extragalactic Database (NED). At the 
adopted distance, the image scale is 1 pixel — 0.05" — 2.2 
parsecs. 



2 GOALS OF THE STUDY 

The most prominent focus of the present discussion is the 
characteristic scale size of the GCs (the effective or half- 
light radius, denoted rh)- We refer to the distribution func- 
tion of their sizes (effective radii) as the SDF. At the 9- 
Mpc distance of M104, a GC with a typical size ~ 3 
pc subtends an angular diameter of about 0.14", which is 
similar to the ACS/WFC point-spre ad funct i on wid th of 
FWHM 0.1". In the sense defined bv [Harris I (|2009ah . the 
clusters are therefore in the partially resolved regime: their 
individual profiles are easily distinguished from stars, and 
their effective radii are clearly measurable. Extensive tests 
of GC profile-fitting for galaxies ranging out to distances 
of ~ 50 Mpc show that GC effective radii can be accu- 
rately measured as long as the stellar PSF is precisely known 
and the cluster size rn is no smaller than about 10% of 
the PSF fwhm (iKundu fc Whitmore I Il99^: iLarsen I Il999l : 



[Jordan et ai~l l2005l : iGeorgiev et al. I l2009l : Harris I l2009al ) 
As will be seen below, the clusters in M104 are far above 
this resolution limit. 

We use our new measurements to establish the SDF 
more accurately, and to define the "Fundamental Plane" 
of GC structural parameters in comparison with other 
galaxies. The half-mass radius rh,m is an interesting 
and important quantity for GC dynamics because it re- 
mains rel atively constant over many internal relaxation 
times fe.g.lSpitzer fc Thuan lll972l:lAarseth fc Heggie Ill998l : 



iBaumgardt et al. II2002I : iTrenti et al. II2007 I) and thus repre- 
sents an intrinsic scale size that was established closer to 
its formation time. As we will discuss in Section 4 below, 
the physical origin of the SDF is likely to be sensitively de- 
pendent on a combination of the star formation efficiency 
(SEE) and the original scale size {r{init)) of the protoclus- 
ters. The fact that the SDF has a near-universal form among 
globular cluster systems can then be turned around to place 
surprisingly close limits on the SEE and r{init). 

The EST cameras, beginning with WFPC2, have pro- 
vided the key instrumental tools that have allowed us to 
resolve GCs in galaxies well beyond the Local Group and 
to measure their scale sizes acc urately. Pioneering work of 



this t ype began a decade ago l|Kundu fc Whitmore 
I2OOII : IXundu et anil999l : iLarsen lll999l : lLarsen e~. 



11996 . 



I2OOI ) 



and the changeover from WFPC2 to the ACS camera soon 
led to faster gains and large databases of cluster sizes 
(e.g. [Jordan et al"~l l2005l . l2009l : [Harris I [2009ah . These sur- 
veys of resolved GC populations have started to reveal 
that rh also depends weakly on three other cluster prop- 
erties. First of these is a g radual increas e in mean 
with galactocentric distance |Hodge I 19621: iMateo 1 1 19871: 
van den Bergh al~l 199ll : Jordan et al. 20051: Spitler et al. I 



2006': Cantiello et al."' '2007*: ' Gomez fc Woodlev I l2007l : 
Harris 2 009a : .Madrid et al. 200^. This trend is generally 
interpreted as the result of the weaker extern al tidal field of 



the host galaxy with increeising distance (s ee von Hoerner 



19571: lKing||l962l: Ivan den Bergh al. Ill99ll: iMurrav fc Lin 



I992I : [Harris fc Pudritz lll994l : lBaumgardt fc Makino II2003I ). 
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raising the intriguing possibility that the GC sizes can pro- 
vide a tracer of the galaxy potential well independently of 
other approaches. 

Secondly, in most galaxies a detectable difference in 
mean size with cluster metallicity has been found, with 
the red GCs being typically 15-20% smaller than the 



blue ones (iKundu fc Whitmore I 



Larsen et al. 



20011: [Jordan et al. 



19981: 
2005| : 



Kundu et aLlllQQgl: 



Gomez fc Woodlevi 



20071 : [Harris I [2b09al . and Paper I). This second trend has 
been variously interpreted as (a) a geometric-projection 
byproduct of the fact that the metal-richer clusters usu- 
ally lie in a more centrally concentrated spatial distribu- 
tion than the meta l-poor ones, coupled wit h the growth in 
mean ru with 7?gc (|Larsen fc Brodie1l2003l ): (b) a result of 
metallicity-dependent stellar evolution ti mes and dyn ami- 
cal mass segregation within the clusters l| Jordan II2004I ): or 
(c) a residual of the cluster formation epoch, in which the 
gas within a metal-richer protocluster could cool and col- 
lapse further t han in a metal-poor one before forming stars 
(|Harrisll2009al '). Combinations of these effects are of course 
not ruled out. 

A third empirically established trend is the correla- 
tion of GC scale size with cluster luminosity or mass. Even 
though there is obvious cluster-to-cluster scatter in intrinsic 
sizes, their median scale size differs little with absolute mag- 
nitude from a global value of 2.5 — 3 pc (e.g. [Jordan et al. I 
I2OO5I : iBarmbv et al~ll2007l ) . It is only at the high-luminosity 
end (~ 5 X lO'^ Solar masses and above) that a trend towards 
somewhat larger rn seems to emerge. This mass/radius re- 
lation may, possibly, bridge the sequence of normal GCs 
to other kinds of compact stellar systems such as UCDs 
(Ultra-Compact Dwarfs) that have scale sizes of typically 
10 - 30 PC jHasegan et al. I l2005l: iKissler-Patig et al 



20061: i Barmbv et al. II2007I : lEvstigneeva et al. II2008I : [Harris 
2009a|). These somewhat more massive compact systems 
(variously called UCDs, Dwarf-Globular Transition Ob- 
jects (DGTOs), or Intermediate-Mass Objects (IMOs); see 
[Forbes et al. [ [20081') obey a different scaling rule ~ M°- ® 
or higher f see [Forbes et al. [[20081 : [Evstigneeva et al. |[2008t ). 
Globular clusters luminous enough to lie in this intriguing 
transition region of 10^ — lO^M© are so rare that it is not 
yet clear whether there is a "universal" mass/radius relation 
for them, or what its detailed form is. 

For all the reasons outlined above, high-quality mea- 
surements of GC scale sizes and structural parameters are 
of great value. Because of its proximity and large GC popu- 
lation, M104 provides one of the few best opportunities for 
this kind of measurement outside the Local Group. 



2.1 The Data and the Profile Models 

Preprocessing of the raw image data is described in Paper I. 
The mosaics of six pointings in each filter were constructed 
at STScI and are publicly available. Paper I also lays out in 
detail the selection of the many hundreds of individual GCs 
in M104 field and rejection of contaminating non-GCs scat- 
tered across the field. The contaminants are mostly small, 
faint background galaxies, but also include any starlike ob- 
jects (rh ~ 0) according to the profile fits. Selection also in- 
cluded the use of color, magnitude, and the distribution of 
objects in a similar control field. There are some 658 identi- 
fied high-probability GCs brighter than R = 24(Mfl ~ -5.8) 
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Figure 1. Location of the target objects within M104 mosaic 
field, which has dimensions 600" X 400". The 659 GC candidates 
are shown as open circles, while the 179 stars used for the PSFs 
are marked as small dots. The large star at bottom shows the 
location of the UCD l|Hau et al. Il2009^ . 



plus one UCD in the final selected list, and we start the 

present paper with that list in hancL 

In Paper I, the ISHAPE code of [Larsen [ (|l999[ ) was used 
to determine the GC effective radii through convolution of 
pXing (1962) model profiles with the stellar PSF. In each 
of the six original ACS/WFC fields a single mean PSF was 
defined, and a profile with a central concentration parameter 
c — Vt jvc = 30 was adopted. For each individual cluster the 
intrinsic profile width th (in practice, the fwhm) is varied 
until a best fit between the convolved profile and the real 
GC is obtained. 

The ISHAPE code is a well tested and fast method for 
obtaining useful distributions for large s amples of par- 
tially or marg i nally r esolved objects (see, e . g. [Larsen |[l99£ ; 



Larsen et ahT [200lj : [Ceorgiev et al. I [20091 : 



Harris I 12009a 



Madrid et alTlbOOgt ). However, it does not easily handle situ- 



ations where the PSF may depend on position on the image, 
and at a more fundamental level it uses only analytic fitting 
functi ons to describe the GC profiles. Although the [King I 
1 19621 ) function was constructed to match real GCs in the 
first place, it does not capture the full range of actual profile 
shapes for GCs, UCDs, and related objects. 

For the presen t work , we use the profile-fitting code of 
[McLaughlin et al~l (|2008l ) . This code has been applied pre- 
viously to GCs in the Milky Way, the Magellanic Clouds , 
M31, and NGC 512 8 iMcLaughlin fc van der Marel I [20051 : 
[Barmbv etall [20071 : [McLaughlin et al. [[20081 ). For the in- 
trinsic GC pro file, this algorit hm can employ either the 
[King [ l|l966l i or [Wilson I l| 19751 ') profiles, both of which are 
true dynamical models constructed from simple but real- 
istic assumptions about the energy distribution function 
of the stars within the cluster. In the discussion below. 
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we refer to these models as "K66" and "W75" Although 
these non-analytic models require a much larger investment 
in computing time, their best-fit solutions lead immedi- 
ately to physically relevant quantities such as relaxation 
times, mean densitie s, binding energies, and fundamental- 
plane parameters ([McLaughlin &: van der Marel ] I2OO5I : 
iMcLaughlin et al. IbOOsf T 

A roughly similar code (KINGPHO T) has bee n used 
for GCs in the V irgo galaxies ( ,Jordan et al. 1 12005| . |2009| : 
iPeng et al. 1 12009|) . Although the King model is quite fa- 
miliar in the literature and describes most real GCs sat- 
isfactorily, the Wilson model has some distinct advan- 
tages in its ability to handle real clusters that have sym- 
metric but radially extended light profiles which (within 
the confines of the King model) would nominally be 
thought to h ave "extra-tidal light". As i s extensively 
discussed by McLaughlin fc van der Marel I Hooi) and 
IMcLaughlin et al. I ( 20081 ) . the kev difference in the Wilson 
model is its smoother ramp-down of the stellar energy distri- 
bution function towards higher energies, which can handle 
relatively larger numbers of stars near the escape energy 
and a more extended envelope with a nominally larger tidal 
radius. 

It is important to note that for small to moderate 
radii within the cluster (that is, out to a few r^), both 
the King and Wilson models give very similar profile shapes 
and do equally well at fitting most real clusters. These in- 
ner regions are populated by the low-energy stars far in- 
side the tidal (escape) radius where the difference in the en- 
ergy distribution function between the two models is small. 
For the clusters in the nearby Local Group galaxies where 
the profiles can be thoroughly sampled over large radii, 
th e Wilson model often performs better, but (as is shown 
in [ McLaughlin fc van der Marel 1 12OO5I : iMcLaughhn et al. I 
bood ) there is no blanket ruling: decisions are made on a 
cluster-by-cluster basis. 

The Sombrero ACS mosaic field does, however, offer a 
particularly vivid case where there is a strong difference be- 
tween these two models. There is a single clearly identifiable 
luminous UCD belonging to the galaxy, whose profile has 
an extended outer envelope which is accurately fit by the 
Wilson model but where the King model comp letely fails. 
We ha ve discussed this special object in detail in lHau et al. I 
I2OO9I). In the present study, we apply both of these two dy- 
namical model profiles to the complete sample of GCs in the 
galaxy and explicitly compare their quality of fit. 



2.2 Measurement of Cluster Profiles 

The profile-fitting c ode i s the same one described fully in 
IMcLaughlin et al. I (|2008l ). For each object, a stellar PSF 
(either an analytic function or an empirical numeric profile) 
is convolved with the K66 or W75 model, and the model 
parameters are adjusted until the convolution fits the ob- 
served surface brightness profile of the cluster. We carried 
out all the measurements on the complete wide-field mosaic 
for the entire field, rather than the six individual fields as 



^ The code can also use any of three analytic fitting functions in- 
cludin g the Sersic profile, a power-law profile, and the equally fa- 
miliar lKingl lll962l) analytic function, which we refer to as "K62" . 



in Paper I. The measurements were done from smoothed 
(median-filtered) images in each filter with the large-scale 
light gradients removed. A 31 x 31 px box filter was used for 
the median filtering. We then used stsdas/ ellipse to measure 
the radial profile for each GC in our list, out to a radius de- 
termined by the point past which the GC surface brightness 
fell well below the sky noise; this limit ranged from > 0.5" 
for the most luminous objects to 0.1" for the very faintest 
ones. 

To define the PSF profile that needs to be paired up 
with each GC, we chose here to select the neares t star to 
each individual GC, as we did in lHau et al. ] l|2009l ). In this 
way we circumvent trying to model the small but potentially 
complex variation in PSF properties across the complete mo- 
saic of six ACS fields. The stars are, in turn, empirically de- 
fined as those objects in our Paper I study that turned out to 
have negligible intrinsic widths (r^ = 0). We drew the PSF 
stars from a candidate list of 179 stars with sufficiently high 
S/N, so the average projected separation between any one 
GC and its accompanying PSF star is ~ 20". The profiles 
of the stars were also measured through stdas/ellipse with 
the same parameters, and each GC/PSF pair was then run 
through the fitting code to find the best-fit solutions. The 
locations of the 658 target GCs and one UCD measured pre- 
viously, along with the 179 candidate PSF stars, are shown 
in Figure [1] As is evident from the figure, we deliberately 
avoid the parts of the galaxy projected on or near the plane 
of the disk, where the profile fits would be badly compro- 
mised by complex background light gradients and crowding. 

Lastly, we re-derive integrated magnitudes B,V,R for 
all the GCs with aperture photometry individually adjusted 
for cluster profile width. T he approach we use is similar in 
principle to lHarrisI (|2009al ) , where in this case we start with 
a fixed-aperture magnitude measured through 5 px radius 
(0.05") and then apply a magnitude correction to "large" 
radius (20 px) Am = ms — m2o that depends in turn on 
the cluster's scale size r^. By direct comparison of the aper- 
ture magnitudes from ISHAPE through these two apertures 
(Figure (21 , we find that Am increases nearly linearly with 
log (rh), 

Am = 0.12 + 0.411ogrh 

and does not differ significantly with bandpass. Thus the 
aperture corrections h ave no signific a nt influence on the 
cluster colors (se e also iHarris 1 l2009al : iJordan et al. I I2OO9I : 
iPeng et al.]|2009l ). 

The output quantities from the proflle fitting code that 
are of primary interest here are the effective radius 
and its uncertainty a-{rh); the central potential parame- 
ter Wo or equivalently the concentration c and its uncer- 
tainty (c = log{rt /rc) where rt,rc are the tidal and core 
radii); and a goodness of fit x^- As all other authors have 
found who h ave worked on GC size measurements for distant 
galaxies (e.g.lKundu fc Whitmore 'l998". '200l"; 'Larson Il999l : 
[Jordan et al. ll2005l : [Georgiev et al. .,2009 : .Harris .200931 ). we 



find that the central-concentration parameter Wo or c is not 
very precisely determined from the best-fit solutions. The 
reason for this is simply that relatively small numbers of 
pixels are available for the code to use for the profile fit 
relative to the PSF size. In turn, some of the other struc- 
tural quantities calculated from c and r^ (such as the cluster 
core radius Vc, which is far smaller than the resolution limit 
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Table 1. Photometry and Structural Parameters for Globular Clusters in M104 



ID 


RA 


DEC 


D 


V 


B-V 


B - R 


rh(pc) 


±(pc) 


c 


± 


1 


12.669456 


-11.639086 


2.777 


18.74 


0.77 


1.26 


3.62 


0.18 


1.73 


0.09 


2 


12.669275 


-11.642228 


2.703 


18.87 


0.83 


1.34 


5.05 


0.06 


1.90 


0.01 


3 


12.666692 


-11.602536 


1.243 


18.83 


0.76 


1.27 


3.72 


0.46 


2.37 


0.71 


4 


12.670719 


-11.641109 


3.874 


18.83 


0.84 


1.36 


4.70 


0.36 


1.94 


0.22 


5 


12.669567 


-11.610695 


2.804 


18.99 


0.80 


1.31 


3.57 


0.10 


1.98 


0.06 




(ISHAPE) (parsecs) 

Figure 2. Empirical aperture corrections to the cluster magni- 
tudes. The mean line through each set of points is given by the 
equation in the text. 



imposed by the PSF) will not be reliable. Empirically, we 
can gauge how reliable our c— values are by comparing the 
output values from the three different filters for the same ob- 
ject. These comparisons show that the internal consistency 
differs widely from one cluster to another, but has a median 
uncertainty near a{c) = ±0.4. 

Fortunately, by far the most robust quantity in the so- 
lutions is the effective radius itself. As is discussed further 
elsewhere (see the references cited above), is relatively in- 
sensitive to changes in c in either direction because, for GCs 
at distances in the ~ 5 — 50 Mpc range, it sits near what is 
essentially a "best" point in the profile: it is neither buried 
in the unresolved core of the cluster, nor lost in the noise 
blanketing its faint outermost envelope. As long as the PSF 
profile is accurately known, any error in c is compensated 
by adjustments in the total profile shape to give the well 
determined radius enclosing half the light. Fortunately, the 
M104 clusters sit well above the empirical resolution limit 
of rh ~ 0.1 f'whm{PSF) below which any leverage on the 
structural parameters is lost. 

A direct comparison of the K66 and W75 models is 
shown in Figure |3l Here, as in iMcLaughlin et al. I l|2008l ) we 
use a normalized ratio defined as (xk — Xw)/{Xj< + Xw) 




R magnitude R ^ (arcmin) 



Figure 3. Left panel: Normalized ratio plotted versus R mag- 
nitude, for the three bandpasses BVR. In this plot, postive val- 
ues indicate that the Wilson (1975) model is preferred to the 
King (1966) model, while values near zero indicate that both are 
equally good. Right panel: Normalized ratio plotted against 
projected galactoccntric distance. 



where Xk,Xw are the values from the King and Wilson 
fits respectively. For most of the clusters, the values of this 
ratio scatter near zero, indicating that both models do about 
equally well at matching the real profiles. Overall, there is 
a slight preponderance for the W75 models to match better 
(positive values in the graphs) but no major differences or 
obvious trends with magnitude are evident. By plotting the 
same goodness-of-fit ratio against projected galactocentric 
distance Rgc, we test for any trends versus background sky 
noise (the large bulge of M104 becomes considerably brighter 
at small Rgc), but we find no systematic trends there either. 

Final successful profile fits and radius measurements for 
652 clusters were obtained. Table [1] lists the results giving 
in successive columns a running ID number; right ascension 
and declination; projected galactocentric distance in arcmin- 
utes; V,B — V, and B — R; the mean rh over all three filters 
and its uncertainty; and the mean King central concentra- 
tion parameter c and its internal uncertainty. The complete 
version of the table is available in the electronic edition. In 
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Figure 4. Some sample profile fits with the King (1966) model 
for six of the clusters, ranging from the brightest in the list to the 
faintest (top to bottom). Dots are the measured surface bright- 
nesses from stsdas/ ellipse, while solid lines are the PSF-convolved 
K66 model solutions. The dashed line at bottom shows one of the 
PSF profiles. 



Figure 5. Left panels: Comparison between the cluster sizes mea- 
sured in the three different filters, as measured from the Wilson 
1975 model (upper panel is B vs. R, lower panel is V vs. R\ all 
scales are in parsec units). Right panels: Same internal compari- 
son for the measurements with the King (1966) model. 



Figure |4j we show several solutions to individual clusters 
with the K66 model and drawn from the _R— band filter. 

2.3 Measurement Uncertainties and Error Budget 

The measurement uncertainties on the rh values were eval- 
uated with a series of tests. First of these was the internal 
consistency in the size measurement among the three filters. 
This comparison is shown in Figure [5] for both the K66 and 
W75 models. In principle, the three filters give three inde- 
pendent measurements of the same quantity for the same 
cluster (and for exactly the same PSF star), with random 
scatter due only to the internal uncertainty of the fit. We 
find good systematic agreement amongst the filters; all three 
correlations scatter closely along the 1:1 lines, so we make 
no systematic corrections between filters. However, the K66 
model fits (right panels in the figure) show slightly smaller 
scatter, fewer outliers, and thus higher internal consistency 
than the W75 model. 

In the end, we find little to choose between these two 
competing models for most individual clusters. Because of 
its slightly superior internal consistency (Figure[5ll , we adopt 
our results from the K66 model for the analysis and discus- 
sion following in the later sections. The K66 reductions also 
allow easier comparison with characteristic-size GC data 
from other galaxies, which we use in the later sections. As is 
discussed above, the W75 model becomes most effective for 
GCs that happen to have extended, low-surface-brightness 
envelopes, which would become noticeable only beyond our 
radial measurement limit of ~ 0.5" ~ 20 pc, and even then 
only for the most luminous clusters whose envelopes would 



be detectable above the sky noise. Higher-S/N data than we 
have at present will be needed to trace these outer parts 
for most clusters (with the notable exception of the UCD 
discussed above). 

To obtain a final set of values, we took an un- 
weighted average of the measurements in all three filters 
and calculated the uncertainty of the mean as equal to 
'^{rh) = {E('"M - {Th)f/n{n - l)Y^^ . The distribution of 
these uncertainties is shown in Figure[6l along with their de- 
pendence on R magnitude. As expected, the average o{rh) 
increases for fainter objects due to lower S/N and the in- 
creased relative effect of background noise. The median un- 
certainty over the entire dataset is ±0.20 pc. 

Next, in Figure [T] we show the difference in the aver- 
age cluster size Ar^ = (rh(K66) — rh(W75)) as functions 
of cluster size and brightness. The median Ar^ is 0.02 pc, 
indicating no important systematic difference between the 
two models over the entire range of the data. There is a 
slight tendency for Ar^ to vary nonlinearly with either size 
or brightness. However, these trends fall within the internal 
scatter, and the most obvious interpretation is simply that 
at some level, we reach the irreducible "floor" where the 
fundamentally different assumptions built into the different 
models can lead to slight differences in the best-fit struc- 
tural parameters. These parameters are unavoidably model- 
defined, and at this level the question about which solution 
represents the "true" cluster size becomes moot. In general, 
we take this graph as useful primarily for estimating the 
internal uncertainty of the fit because of the model assump- 
tions alone. This fitting uncertainty per cluster (which we 
adopt as the rms scatter in the graph) is then ±0.21 parsec. 

A final internal check of our reductions is to gauge the 
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Figure 6. Upper panel: Histogram of uncertainty (T(r^) in cluster 
effective radius returned by the model fit, in parsec units. Lower 
panel: Fitting uncertainty versus cluster brightness. The median 
uncertainty of 0.20 pc is shown as the dashed line. 



Figure 7. Left panel: Difference between the cluster sizes Ar^ = 
(r^ (King)— (Wilson) plotted versus size. Both scales are in pixel 
units (1 px = 0.05"). Right panel: Ar^ versus cluster brightness. 



uncertainty in size measurement due to the PSF itself. Our 
"nearest star" approach to defining the PSF for each target 
cluster ensures that the results will not be affected by large- 
scale trends in PSF size across the mosaic, but any one PSF 
star will have lower S/N than the average of many of them 
across the field and thus slightly higher random uncertainty. 
To evaluate this level of uncertainty, we ran the model fits 
a second time, now using the second-nearest star for each 
object. Direct comparison of the two reductions is shown 
in Figure [S] Again, no systematic difference appears bigger 
than 0.05 pc, but the rms scatter is ±0.30 pc. 

In summary, all three of the possible sources of fit- 
ting uncertainties (filter-to-filter consistency, choice of pro- 
file model, choice of PSF) contribute to the net uncertainty 
to about the same level. Adding the three in quadrature, we 
then estimate that a global-average uncertainty in the clus- 
ter sizes, due strictly to the internal measurement process, 
is ±0.41 parsec. As will be seen below, this is equivalent to 
about 16 percent of the median cluster size. 

It is worth noting here that the mean uncertainty of 
±0.4 pc also allows us to resolve an important feature of the 
GC size distribution, namely its lower limit. As will be seen 
in Section 3, the SDF starts increasing sharply near Vh ~ 1 
pc; clusters smaller than this are almost nonexistent in M104 
or in other galaxies. For the Virgo galaxies at d = 16.5 Mpc, 
almost twice as far away as M104, this lower edge to the 
SDF is not well determined I Jordan et al. pOOSi) except with 
extremely high-S/N data ^Madrid et al. 200^). For the gE 
galaxies at d > 30 Mpc studied by Harris (2009af), the edge 
falls below the limits of HST resolution and only the upper 
half of the SDF can be clearly measured. The relevance of 
this feature to understanding the formation of GCs and the 
origin of the SDF will be discussed in Section 4 below. 
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Figure 8. Difference between the measured cluster sizes for two 
different assumptions about the PSF star. One run (denoted here 
PSFl) uses the nearest moderately bright star to each cluster as 
the fiducial PSF for that cluster; the second run (PSF2) uses the 
second- nearest star. 
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2.4 Comparisons with Paper I 



Lastly, we compare the results from our new model profile 
fits with those done in Paper I, which used the ISHAPE code 
l|Larsen Ill999l ) with the analytic K62 model profile. For this 
test, w e carried out a separate run of the McLaugh lin et al. I 
l|2008ll code used above, but now using the K62 model in 
order to rule out any differences due only to the adopted 
model and focus on the two codes themselves. 

In the first panel of Figure (5] we show the direct com- 
parison of this run with our K66 fits. Though there is overall 
close agreement, the median difference is 0.19 pc in the sense 
that the K66 models tend to measure the clusters slightly 
(8 percent) larger. The rms scatter is ±0.32 pc, similar to 
the internal comparisons already described above. 

The second panel of Figure [9] shows the correlation of 
our new K62 fits with those from ISHAPE and Paper I. For 
small objects {rn < 2.5 pc, or effective radii less than about 
1 pixel on the images) the two methods agree systematically 
quite well. However, for larger objects (> 2.5 pc) there is an 
offset between the two codes in the sense that ISHAPE mea- 
sures them smaller by about 0.3 pc than does the McLaugh- 
lin et al. code. One difference between the two runs is that 
the ISHAPE reductions assumed a fixed central concentra- 
tion of c = log(rt/rc) = 1.5 for all clusters, whereas the 
McLaughlin code solves for c (equivalently, the central po- 
tential Wo) as a free parameter. However, this should mainly 
introduce cluster-to-cluster scatter since c ~ 1.5 is a reason- 
able average for real GCs. A possible cause for a systematic 
offset may be (as is described in Paper I) that the ISHAPE 
sizes in the three filters were all normalized to the V^— band 
data and then averaged, but it is not clear why a magnitude- 
dependent offset should appear. These two factors aside, the 
remaining differences are presumably due to the details of 
the two codes, and we use this test only as a rough con- 
sistency check. We conclude that at the 10-20% level, both 
codes and all three models return similar intrinsic cluster 
sizes. As will be discussed in the next section, most trends 
found in Paper I, such as cluster size versus magnitude or 
galactocentric radius, fall into the same patterns here. 



3 THE SIZE DISTRIBUTION 

The color/magnitude diagrams in {R, B — V) and [R, B — R) 
are shown in Figure [TO] The classic bimodal division into the 
blue, metal-poor and red, metal-rich sequences is easily vis- 
ible, with a (more or less arbitrary) division at (B — R) ~ 
1.30. A closer analysis of these bimodal sequences is pre- 
sented in Section 5 below; first, we take a closer look at 
the distribution of GC sizes as a whole, and their correla- 
tion with external properties including luminosity, metallic- 
ity, and projected galactocentric distance. 
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Figure 9. Left panel: Cluster effective radius from the King 
(1966) profile fits, plotted versus our average r^j values determined 
from the simpler King (1962) analytic model. Right panel: Cluster 
effective radius as determined from ISHAPE in Paper I, plotted 
versus the King (1962) model fits determined in the present paper. 
A 1:1 correlation line is drawn in for each graph. 
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3.1 The Overall Scale Size Distribution 

The overall distribution of the GC scale sizes in our dataset 
is shown in histogram form in Figure [TT] The distribution is 
modestly skewed to larger rn, with a median at 2.44 ± 0.04 
pc and an intrinsic dispersion of ±0.85 pc. 

To calculate the dispersion for this and for other dis- 
tributions used in this paper, we adopt here the Median 



Figure 10. Magnitude R versus color index {B — V) or [B — 
R) for 659 globular clusters in M104. The sin gle known UCD i s 
shown at top as the large starred symbol fsee lHau et al. Il2009h . 
All BVR magnitudes are those measured through 5-px aperture 
individually corrected to 20-px radius for cluster size as described 
in the text and Figure [2] 
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Figure 11. Histograms of globular cluster scale sizes for three 
galaxies. Our present results for M104 are shown as the shaded 
histogram; the Milky Way GCs are in shown as the dashed red 
histogra m, and the compilatio n of old GCs from nearby dwarf 
galaxies ijCeorgiev et al. 1 12009 ) as the dotted blue histogram. The 
Milky Way and dwarf samples have been normalized to the same 
total population as in M104. The errorbar at upper right is the 
mean internal uncertainty of ±0.41 pc (see text) for our measure- 
ments. 

Absolute Deviation (MAD), a robust estimator of the in- 
trinsic scatter of a data sample. It is useful in cases where 
the distribution is asymmetric and even for cases where the 
conv entional standard de viation may be formally undefined 
(e.g. iHoadin et aTlllQSd ). For a dataset {xi} with median 
X the MAD is defined as 

MAD = median{|a::j - x\} (1) 

and the sample dispersion is then estimated as 

a = 1.483 X MAD. (2) 

For a Gaussian distribution, this formula exactly gives the 
usual standard deviation. 

A basic point of immediate interest is to compare the 
GC size distribution with the one for the "baseline" Milky 
Way system. However, a minor spatial bias must be kept 
in mind. The Sombrero GCs in our list do not cover its en- 
tire halo, whereas (in the Milky Way at least) there is a 
well known trend for (r^,) to increase systematically with 
Rgc (|van den Bergh al. Ill99l1 ). The GCs in our data have 
projected galactocentric distances ranging from Rgc ~ 0.9 
kpc out to 15 kpc, with reasonably complete radial coverage 
to 8.5 kpc. To extract a Milky Way sample that will more 
closely mimic the M104 data, we take the 114 known Milky 
Way clusters wi th measured half-light radii obtained from 
King- model fits (jHarris Ill996l ) and with projected Galacto- 
centric distances Rgc = {Y^ + Z^)'^^^ < 12 kpc. Here {Y, Z) 
are the distance components projected on the sky paral- 



lel to and perpendicular to the Galactic planelf] In Figure 
111! the distribution for the 114 selected Milky Way clus- 
ters is shown as the dashed red histogram. The median is 
at rh{MW) = 2.61 ± 0.10 pc, and to first order, the two 
distributions are strikingly similar. The intrinsic spread of 
cluster sizes around the peak is characteristically ~ ±1 pc 
rms, and very few clusters exist with physical scale sizes 
smaller than about 1.0 pc. The slightly broader width of the 
central peak for M104 is likely to be due in part to the in- 
strumental broadening of ±0.4 pc. (However, it should also 
be noted that the Milky Way data come from a somewhat 
heterogeneous collection of starcount and surface brightness 
data obtained in many programs for clusters at very differ- 
ent distances from the Sun). 

Another dataset making an interesting comparison is 
the recently measured set o f old G Cs in several dwarf galax- 
ies, from Georgicv et alT] (|2009l ). On average these small 
galaxies are at similar distances to M104, and their profiles 
were meas ured from HST/ACS iin ages with ISHAPE±K62 
profile fits. iGeorgiev etlTl (|200i) give data for 83 objects 
regarded to be classically "old", luminous GCs. Of these, 
42 clusters come from dirr galaxies, 31 from dSph and dE 
systems, and 12 from Sm galaxies. The size distribution of 
these is plotted in Figure [TT] as the blue dotted histogram. 
The median of this dwarf-galaxy sample is at = 2.87±0.19 
pc. 

Although the medians or the peak points of these three 
samples are not strongly different, a potentially more impor- 
tant test is the total shape of the whole distribution. A stan- 
dard Kolmogorov-Smirnov two-sample test shows that the 
M104 and Milky Way GCs are significantly different at the 
97% level, whereas the M104 and dwarf samples are differ- 
ent at more than 99% confidence (the Milky Way and dwarf 
samples are not significantly different from each other) The 
key difference among these histograms is the relative num- 
ber of "extended" clusters with > 5 pc, i.e. the ones more 
than twice the median size. Our initial selection of GC can- 
didates in M104 did not rely on object scale size except for 
rejection of small, starlike objects (see Paper I), so the sam- 
ple should not be biased against GCs in the range ~ 5 — 10 
pc or even larger. Nevertheless, even a few more objects 
added to the high-end tail of the distribution would notice- 
ably reduce the statistical difference between the Milky Way 
and M104, so for the present, we regard these comparisons 
as only indicative. Perhaps a more important conclusion is 
that if we use only the clusters smaller than 5 pc, there are 
no significant differences among these three samples. 

Extended clusters would clearly find it easier to sur- 
vive in the gentler tidal environment of smaller galaxies, or 
in the outer-halo regions of large galaxies. If we take the 
comparisons in Figure [9] at face value, M104 - a massive, 

^ The third component X is directed along the axis from the 
Sun to the Galactic center, and contains most of the random 
erro rs in the distance meas urements to the individual clusters 
fsee iRacine fc Harrisi Il989l) . By projecting them onto the YZ 
plane we therefore get the closest to seeing the system as if it 
were an external galaxy free of position-measurement bias. 
^ The difference between the M104 and dwarf samples would be 
even stronger if, as is hinted by Figure (9] the ISHAPE fits need 
to be corrected to slightly larger r/i to put them onto the same 
internal scale as the present code. 
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bulge-dominated Sa galaxy - has very few such extended 
GCs compared with the other two samples. Possible inter- 
pretations that immediately suggest themselves are either 
that M104 did not acquire most of its globular cluster pop- 
ulation by accretion of small satellites; or that any extended 
clusters that might have been accreted this way have already 
been tidally destroyed. A survey of the outer parts of M104's 
halo would provide much clearer evidence to discuss this ar- 
gument further. 



3.2 Correlations with Metallicity 

Previous surveys of GCs in nearby galaxies have shown that 
the blue, metal-poor clusters are consistently lar ger on av- 
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[Jordan et al. i .2005: Harris 2009al ) they are highly signifi- 
cant in a statistical sense, a nd the y are found at all galac- 
tocentric distances. [Harris I (|2009al ) determines a mean size 
difference of (17±2)% (blue minus red) from a sample of sev- 
eral thousand clusters in six supergiant ellipticals (Brightest 
Cluster Galaxies or BCGs). 

Our observed correlation is shown in the upper panel of 
Figure [121 The red clusters do lie at lower sizes on average, 
but an unweighted linear fit to the complete dataset gives a 
slope Arh/A{B - R) = -0.245 ± 0.191, which is not highly 
significant. (See also Paper 1 for a similar diagram.) For all 
blue clusters (defined here as those with {B — R) ^ 1.30), 
the median size is rhiplue) = 2.47±0.04 pc, while the redder 
clusters (B — R) > 1.30) have a median rhired) = 2.33±0.05 
pc. The metal-poor GCs are thus 6% larger, and the differ- 
ence Ar^ = 0.14 ±0.064 pc is significant at the 2.2a level, a 
tantalizing but not conclusive offset. In Paper I, the differ- 
ence between median sizes was found to be Arji = 0.25±0.06 
pc, similar to the present value. Thus we find the same effect 
seen in other systems, but at a smaller amplitude. 

To date, the strongest claims for metallicity-related dif- 
ferences in size have relied on the large GC samples from 
elliptical galaxies. The samples from GCs in disk galax- 
ies are much smaller and inev itably have lowe r stati stical 
significance. Notably, however, [Cantiello et al. [ l[2007l ) and 
[PeGraaff et al. [ l|2007l ') find that the mean size difference for 
the GCs in the disk galaxies NGC 5866 and 1533 is at the 
level of ~ 0.1 pc or less, quite similar to our results for M104. 

Correlation of mean size with location in the halo has 
been proposed as an important factor in deciding what 
the physical cause for the metallicity /size offset actually is 
l|Lars cn & Brodic 2003; Jordan 2004). If the size difference 
versus metallicity persists at all galactocentric distances, 
then it is less likely to be due to a geometric projection 
effect (the metal-richer clusters are usually found to have 
a more centrally concentrated spatial distribution than the 
metal-poor ones, thus will be more subject to stronger tidal 
stripping. This effect would be much stronger on the inner- 
halo clusters, and not as important for the clusters of both 
types that are found in the outer halo). In M104, as shown 
in Paper I, both red and blue GC subsystems accurately 



For the red GCs, using our present data over the radial re- 
gion Rgc < 3' where we have complete azimuthal coverage, 
we find a = -1.937 ± 0.046,/? = 3.191 ± 0.053; while for 
the blue GCs, a = -1.335 ± 0.179, /? = 2.541 ± 0.205. Thus 
the blue, more metal-poor subsystem is significantly more 
extended. According to the geometric projection hypothe- 
sis, we would then expect a significant mean size difference 
between blue and red in the inner regions but less so in the 
outer halo. 

This version of the size distribution is shown in the 
lower panel of Figure 1121 separately for the blue and red 
clusters. Median values for are plotted as the connected 
large points with errorbars, in 1-kpc bins. The GCs of both 
metallicity groups have scale sizes that increase gradually 
but consistently throughout the halo of the galaxy. The 
slopes of both trends are similar, so if we combine all clus- 
ters to gain statistical weight, we find a rate of increase of 
mean cluster size ~ Rg'^'^^'^''^^ . The effective radius of 
the spheroid light is Reff ~ 1.2' = 3.14 kpc, so the to- 
tal radial range of our data reaches to an outer limit of 
about 4.5-Re//- Not only does the radial increase affect all 
clusters, but the difference between red and blue remains 
similar (and small) at all radii. This large-scale trend may 
therefo re provide evidence aga inst the geometric-projection 
effect (|Larsen fc Brodiel [2003[ ) as being the sole explana- 
tion for the size difference in this galaxy. However, a more 
detailed deprojection model will be needed to test th is con- 
clusion more quantitatively (jLarsen fc Brodie1[2003h . 

The only other data in the literatur e that cover similarly 
large ranges in Rgc are for NGC 5128 ([g omez fc Woodlev I 
[2003), which reach even further to ~ 8Reff but cover a 
sm aller sample; an d for the six supergiant ellipticals studied 
bv [Harrisi l[2009al ). which extend to distances ~ 4.5i?e//. 
Both of these other surveys, however, indicate the same 
steady increase in scale size with Rgc for clusters of both 
types. So do the s maller samples in the two disk galax- 
ies m e ntioned above ([Cantiello et al. |[2007l : [PeGraaff et al. I 
[2007^ . [Harrisi (|2009a[ ') derives a simple power-law scaling for 
gE's of rh{med) = rho{Rgc/ Ref f)'^'^^ , where the zeropoints 
are rno = 2.53 pc (blue) and 2.15 pc (red). These functions 
are shown in Figure [T2T b). Clearly, they strongly resemble 
our current data, bracketing the median points and indicat- 
ing that the M104 GCs follow a very similar trend. 

The total evidence suggests that the dependence of GC 
scale size on metallicity is at least partly intrinsic to the 
clusters, and thus due to a more local cause having to do ei- 
ther w ith their formation or later internal evolution . [ JordaiTI 
([2004h has proposed that it is the result of stellar-evolution 
timescales that depend on metallicity, coupled with many 
internal relaxation times of dynamical evolution that would 
make the metal-rich clusters appear smaller. 

Another possibility is simply that the metal-richer clus- 
ters benefitted from more rapid cooling and contraction 
while they were still gaseo us protocluste rs and had not yet 
formed most of their stars ([Harris 1 [2009a[ ). Still another pos- 
sibility is that all clusters started out with similar sizes dur- 
ing their early protocluster stage, but the star formation 
efficiency (SFE) was a bit higher for more metal-rich gas, 
allowing the cluster to expand less during the gas expulsion 
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Figure 12. Upper panel: Scale sizes of M104 clusters versus {B — 
R) color. A linear fit to the points is shown as the dashed line. 
Lower panel: Scale sizes for the clusters as a function of projected 
galactocentric distance Rgc- Blue-sequence clusters [B — R^ 1.3) 
are plotted as crosses, red-sequence clusters [B — R > 1.3) as open 
circles. Large filled circles with errorbars denote the median sizes 
in 1-kpc bins, separately for the blue clusters (dashed blue line) 
and red clusters (solid red line). The two smooth curves (blue, 
red) show the scaling laws r^ ~ ^^'c^ that have been found for 
globular clusters in supergiant ellipticals jHarris Il2009a^ . 



Figure 13. Upper panel: Sizes of the blue-sequence clusters 
{{B — R) < 1.3) as a function of luminosity M^. The connected 
solid dots with errorbars show the median rf^ in half-magnitude 
bins (see text). Lower panel: Sizes for the red-sequence clusters 
{{B — R) > 1.3). In bot h panels, the m ean relations for GCs in 
six supergiant ellipticals llHarris i boOQal) are superimposed as the 
dashed lines. 



phase (see Section 4 below). These alternatives will be dif- 
ficult to compare quantitatively, but detailed models would 
be of great interest. 



3.3 Size versus Luminosity 

The correlation of cluster scale size with luminosity is the 
observational version of the mass/radius relation. This form 
is displayed in Figure 1131 Here, the individual clusters are 
plotted along with the median in half-magnitude bins of 
luminosity Mr — R — (m — M)r where our adopted ap- 
parent distance modulus is (m — M)_b — 29.90. Over the 
range —7 > Mr > — 10 (corresponding approximately to 
the luminosity range 3 x 10^ to 6 x 10^ Lq) the median size 
remains nearly constant, increasing more rapidly for lumi- 
nosities Mr < —10. The pat t ern for the GCs in the super- 
giant ellipticals from [Harris I (|2009al ). shown as the dashed 
lines in Figure lTS] is the same, although these six galaxies are 
~ 5x more distant than M104 and thus the GCs could not 
be traced to similarly faint luminosities. Here, the (7, B — I) 
photometry used for the giant ellipticals has been converted 
approximately to Mr with the assumption R — I ~ 0.5. 

In Figure [T4l the data for M104 are plotted in compar- 
ison with two other l arge disk galaxies, th e Milky Way and 
M31 (with data from lBarmbv et al. |[2007l '). This form of the 
graph shows more clearly the trend for the most luminous 
GCS (L > 3 X IO^Lq), extending up to the very most lumi- 



nous GCs known near lO^'L©. A more extensive discussion 
of the possible link between massive GCs and the still more 
massive UCDs is made by iHasegan et al. 1 ll2005h and has 
been develope d further in, f or example , 'Kisslcr-Patig et al. I 
(2006); Evstigneeva et al. I l|2008l ): iForbc s ct al. (200i). In 
general, UCDs have scale sizes of ~ 10 pc and above, thus 
sit a bit higher on the (L, r^) relation than do the GCs. The 
single UCD known in the M104 field is more than one mag- 
nitude brighter (Figure llOp than the top end of either the 
blue or red GC sequences. 

A key factor distinguishing a massive GC from a UCD 
or nuclear cluster ma y be the h i gher mass-to-light ra- 
tio fo r UCDs (Micskc et al. I l2008l : iBaumgardt fc Mieskel 
|2008| ). Multiple cluster mergers are another route to form- 
ing supermassive GCs and UCD-like objects and can ex- 
plain the observed uptur n in the L vs. correlation (e.g. 
iKissler-Patig et al. I [20061 ): but it is less clear whether such 
mergers would also produce objects with increased (M/L). 
It remains possible that at least some UCDs may s imply be 
very m assive GCs. This latter view is supported by [Murray I 
(|2009l ). who develops a model for the size/mass/luminosity 
relation using the idea that very massive protoclusters (> 
10^ Mq) will be optically thick to their far-IR radiation, 
making radiation pressure important for energy balance. 
The resulting increase in Jeans mass with cluster mass will 
yield a scaling of M/L with mass that matches the trend for 
UCDs and the most massive GCs reasonably well. 
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Figure 14. Scale sizes of globular clusters in three massive disk 
galaxies: M104 (small crosses), the Milky Way (large filled circles), 
and M31 (triangles). The median size of 2.4 pc is drawn in as the 
dashed line. Note the systematic rise in median rj, for luminosities 
above log [Ly /Lq) > 5.5, as well as the scatter to high at low 
luminosities. 



3.4 Fundamental-Plane Quantities 

Recent work has established the existence of a surprisingly 
narrow "fundamental plan e" (FP) of structural quantities 
for g l obular clusters (e.g. Diorgovski I 19951 : iMcLaughlin I 
I2OOOI : iMcLaughlin fc van der Marel I bOOsf T The so-called 



K— space of three orthonormal quantities (ki, K2, K3) is con- 
ventionally constructed from the central velocity dispersion 
(jpo, effective radius r^, and surface mass density E^. For 
GCs in galaxies such as this one beyond the Local Group, 
direct velocity dispersion measurements (hence mass) are 
rare by comparison with photometry, and in any case require 
knowing the cluster core radii rc to convert the integrated 
velocity dispersion to its core value. The core radii in turn 
have to be estimated from and the model-fitted central 
concentrations c, which are quite uncertain for objects such 
as these where is far smaller than the image resolution. 

Although a full discussion of the FP in all its vari- 
ous guises is therefore not possible, more limited versions 
can still be constructed. A useful quantity representing the 
cluster density is the surface intensity Ih ~ I/tot /(27rr^) 
integrated over the half-light radius. This quantity is well 
enough determined to allow us to compare the M104 GC 
pop ulation directly with the baseline Milky Way system (see 
also lBarmbv et"aL l l2007l . for similar data in M31 and other 
Local Group members). Our measurements give internal un- 
certainties of ±0.1 mag or better in luminosity and ±0.4 pc 
in radius, making (log Ih) uncertain to ±0.13. In Figure 
[T5I Ivh is plotted against both cluster size and total lu- 
minosity. The analogous Milky Way data are taken from 
IMcLaughlin fc van der Marel I (|2005i). The obvious trends 
that Ih is brighter for more luminous or more compact clus- 



Figure 15. Surface intensity in V over the half-light radius, plot- 
ted as a function of cluster scale size (left panel) and total lu- 
minosity (right panel). Small dots give the data for the M104 
clusters, while the l arge red crosses give similar data f or the 
Milky Way clusters l|McLaughlin fc van der Marel 1 12005^ . The 
large square at upper left in each panel with errorbars shows the 
mean internal uncertainty in each datapoint. The straight lines in 
each graph show the nominal proportionality scalings ~ r"'^ 
and Ih ~ L. 



ters are evident from the graphs, but the main conclusions 
to draw are that the M104 GCs match up well with both 
the mean positions of the Milky Way clusters and the in- 
trinsic scatter around the fiducial scaling lines (shown in 
the graph). The clusters at lowest luminosities and surface 
brightnesses do, however, tend to have systematically larger 
scale sizes, dropping them below the fiducial scaling lines 
in both graphs. The papers listed above give more detailed 
discussions of these relations. 

An alternate and conceptually simple formulation of 
the FP is in terms of the binding en ergy of the cluster, 
Eb = f{c)GM^/R (IMcLaughlin I I2OOOI ). where /(c) repre- 
sents the details of the internal mass distribution for a given 
cluster and R is som e ch aracteristic rad i us. As is shown in 
iMcLaughMnI ()2000l ) and iBarmbv et al. I (|2007[ ). for normal 
globular clusters the binding energy varies almost exactly 
as Et ~ with remarkably little relative scatter, a scal- 
ing law that is basically quite different from the Et ^ AI^'^ 
rule characterizing other structures such as giant molecu- 
lar clouds or E galaxies. This form of the FP has been 
used for globular clusters within Local Group galaxies in- 
cluding the Milky W ay, M31, and the Magellanic Clouds 
l|Barmbv et al. II2007I ). As originally defined by McLaughlin, 
the calculation of Ei, requires fairly precise knowledge of the 
King core radius and central concentration, which we do not 
have for M104 and othe r galaxies at large distances. How- 
ever, iMcLaughlinI l|2000l ) also shows that it can be recast in 
terms of the effective (half-light) radius, which is much more 
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accurately known. The major advantage of defining Et this 
way is that it allows us to extend the FP discussion to the 
huge cluster populations in the giant galaxies that lie in the 
near-field region beyond the Local Group. 

An appropriate combination of the equations in 
McLaughlin's paper gives 



= G(-) ^ — 



(3) 



where TZ,£,C are all dimensionless fu nctions of c that ca n 
be calculated from the prescriptions in iMcLaughlin I l|2000l ). 
Happily, the ratio {TZ£ / C ) is nearly constant (see also 
lMcLaughlin1l2000l ). equalling 0.17 ± 0.02 over the range of 
King c— values appropriate for normal GCs, so by using the 
half-light radius we do not need to know the central con- 
centration precisely. A potentially more important caveat is 
that strictly speaking, knowing Ad requires also knowing the 
cluster mass-to-light ratios. In turn, measurement of M/L 
independently of the photometry or colors requires direct 
measurement of the internal velocity dispersions, which are 
not yet available for M104 with the exception of the bright 
UCD (see below). For the present purposes, we simply as- 
sume {M/L)v = 2.0 and then calculate E,, for all the M104 
clusters in our list. This fiducial [M/L) is a compromise 
choice drawn from the recent literature for dynamically mea- 
sured masses of GCs in the Milky Way, M31, and NGC 5128 
jMcLaughlin fc van der Marel I (20051: iReikuba et al. II2007I : 



IStrader et al. I l2009l : iBaumgardt et al. II2OO9I ). These mea- 
surements fall in the typical range {M/L)v — 1 — 3, probably 
with real cluster-to-cluster differences depending on metal- 
licity, luminosity, or galactocentric distance. 

In Figure [161 the resulting correlation of binding energy 
with cluster luminosity is shown. A direct least-squares so- 
lution gives log Bb/ergs = (41.11±0.07) -h (1.92±0.02)logL, 
with a scatter of ±0.16 dex. For the Milky Way GCs 
McLaughlin found slopes in the range 1.8 — 2.2 under var- 
ious assumptions. In summary, we concur that the scaling 
rule Eh ~ provides an excellent first-order description of 
the data. The scatter around the line in Figure 116b , is artifi- 
cially small (only half as large as for the Milky Way sample) 
because it does not account for cluster-to-cluster differences 
in the M/L ratio that must be present (and since E^ varies 
as L^ it is moderately sensitive to changes in M/L). 

The UCD in the M104 field is of special interest as a pos- 
sible "connector" to the GC sequence. Its position is shown 
in Figure [16] at upper right. The lower of the two connected 
points is where we would have located it if we had assumed 
(M/L)v = 2.0 as we did for all the GCs. The upper point 
uses its actual value of {M/L)v = 4. 36 as directly mea sured 
from its internal velocity dispersion (|Hau et al. II2009I ). Re- 
markably, its true position on the graph extends the same 
Eb curve defined by the GCs accurately upward by another 
order of magnitude beyond the top of the GC sequence. This 
result, in ad dition to the othe r characteristics of the UCD 
measured bv lHau et al. I (|2009l ). is consistent with the inter- 
pretation that this luminous, compact system is a massive 
GC. 

The binding energy also provides a sensitive confirma- 
tion test of change s in cluster structure with environment. 
iMcLaughlinI l|2000l ) found that the residuals from Et vs. L 
were a significant function of galactocentric distance in the 
Milky Way (see his Fig. 10) , in the sense that clusters further 
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Figure 16. Upper panel: Binding energy Ej, measured in ergs, 
plotted against cluster V— band luminosity in Solar units Ly/Lq. 
The straight line is an exact scaling law Ef, ^ . The two large 
starred points at upper left show the position of the UCD under 
two different assumptions about its M/L ratio (see text). Lower 
panel: Normalized binding energy E^/L'^ versus projected galac- 
tocentric distance. The fitted straight line has a slope of —0.19 
(see text). 



out in the halo are less tightly bound because of their sys- 
tematically larger radii. Using the results from Figure [T2l and 
Section 3.2, we can predict that for our M104 data, {Eb/L^) 
should vary as ^Rj^"'^'^*""^' . The consistency test is shown 
in Figure [T5b . where we find a net downward trend with a 
fitted slope of (—0.19 ± 0.02) in good agreement with the 
prediction. For the Milky Way, McLaughlin found a slope of 
—0.4 ± 0.1 for the trend of E^/L'^ versus three-dimensional 
Galactocentric distance, which after projection to two di- 
mensions will decrease the slope to a level closer to our re- 
sult. 

Yet another way to represent the trend of cluster struc- 
ture with galactocentric distan ce is as used r e cently by 



ture witn gaiactoce r itric aistan ce is as usea r e cently Dy 
[McLaughlin fc Fall I ([2OO8I ) and [Chandar et aTl (|2007l ). A 
characteristic mean internal mass density for each cluster 
can be calculated from ph = {M/2) x 3/(47r7?|) where 
Rh — (4/3)rh is the three-dimensional half-mass radius. 
As above, we use M/L = 2 to transfer from luminos- 
ity to mass. The results plotted separately for the blue 
(metal-poor) and red (metal-rich) clusters are shown in Fig- 
ure 1171 The best-fit lines through each set of points are 
\ogPh{MQpc-^) = (3.35 ±0.10)- (0.51 ±0.12)log7i9c (blue) 
and logph = (3.33 ± 0.10) - (0.49 ± 0.1S)logRgc (red). The 
scatter around both relations is ±0.55 in log ph- If mean 
GC mass does not vary with location in the halo, as the 
M104 data show, then the scaling ~ Rgc^ (Section 3.2) 
would predict ph ~ RJ^'^^, which m atches what is found in 
the density plot. As is d i scusse d by [ Chandar et al. [ ()2007[ ) 
and [McLaughlin fc Fall [ l|2008l ). the large cluster-to-cluster 
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Figure 17. Calculated mean density plotted as a function 
of galactocentric distance. The metal-poor and metal-rich clus- 
ters are plotted in separate panels, with the mean relationships 
derived in the text shown as dashed lines. 



scatter in density at all galactocentric distances allows the 
mean cluster mass to be nearly independent of Rg^ in the 
presence of density-dependent dynamical evolution times. 



4 ORIGIN OF THE SIZE DISTRIBUTION 

Aside from minor trends with metallicity, galactocentric dis- 
tance, and mass, the characteristic scale sizes and their dis- 
tribution function (the SDF) in all types of galaxies are re- 
markably similar over an impressive range of environments. 
Early indications of the near-universal mean or median GC 
siz es on observational grounds w ere developed a decade ago 
by iKundu fc Whitmore I (|200ll ) from HST/WFPC2 mea- 
surements of GC s izes in 28 elliptical galaxies, and by 
iLarsen et all (|200ll ) for a larger sample of GCs in 17 el- 
li pticals. A major step forward was taken with the work 
of I Jordan et al. I l|2005h . who used their database of ~ 10^ 
GCs in dozens of Virgo galaxies to strongly reinforce the 
view that the full shape of the SDF, as well as its mean or 
median, is a near-universal characteristic of GCs in galaxies 
generally. Furthermore, the SDF has also b een found to be 
similar over a wide range of cluster age (e.g. iBarmbv et al. I 
[2OO6; Schccpmakcr et al. 2007). Its key features are a typi- 
cal median size ~ 2 — 4 pc, a rather sharply defined cutoff 
below 1 pc, and an asymmetric tail extending to larger radii. 
The near-universality of this distribution across all types of 
galaxy environments suggests to us that the conditions local 
to the clusters themselves, during their formation period, are 
an importa nt factor in de t ermin ing the observed size distri- 
bution. As [Jordan et all l|2005l ) remark, "The form of the 
distribution ... should serve as a useful constraint for models 
of GC formation ... any viable picture of star formation in 



clusters should produce an observed size distribution that is 
consistent with the form of [the SDF]". 

A more well known feature of GC systems within 
galaxies in general is their mass (luminosity) distribution 
(GCLF), which is also a near-universal function relatively in- 
sensitive to host galaxy size or environment. The GCLF can 
be modelled as the outcome of a power-law-like initial mass 
function coupled with many Gyr of dynamical evolution, 
which preferentially removes the low-mass and low-density 
clusters and reduces the IMF to the peaked Schechter-like 
form seen today (for only a handful of the dozens of papers 
discussing t he mass distribution function and its secular evo- 
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2009f). Bv contrast, the origin of the linear size distribu- 



tion for GCs has received less attention. But the state 
of development of both observations and models is now 
reaching the point where some paths to understanding it 
are opening. 

As stated above, the GC effective radius remains nearly 
invariant over the normal dynamical evolution of the clus- 
teiQ - much more so than the total mass of the clusters, 
which decreases by factors of 3 or more over a Hubble time 
as it loses stars through the slow processes of tidal strip- 
ping and evaporation. Thus gives us an unusually di- 
rect glimpse of its characteristic size much closer to its ori- 
gin. However, extrapolating the current effective radius all 
the way back to its protocluster epoch is highly unlikely 
to be correct. A crucial stage in its evolution occurs just 
after formation, when it consists of a mixture of gas and 
stars in proportions determined by the star formation effi- 
ciency (SFE). Especially over its first ~ 30 Myr, a cluster 
experiences rapid mass loss due to the energy output from 
its massive stars, including UV radiation, supernova ejec- 
tion, stellar winds, and even external dynamical heating. 
These effects expel the gas and drive an internal expansion 
of the cluster. For the extreme case of low SFE and instan- 
taneous gas loss this stage would lead to rapid dissolution of 
the cluster into th e field. For bound star clusters, however, 
both observations (|Mackev fc Gilmore ^*2003': 'Bastian et al. I 
l2008h and theory (see, e.g. Baumgardt fc Kroupa ,, 20071 . for 



a thorough overview) show that for plausible SFEs, the clus- 
ter expands typically by a factor of 2 — 5 over this cru- 
cial initial stage. An empirical expre ssion for the growth of 
core radius with age for real clusters l|Bastian et al. |[2008l ) is 
rcore = 0.61n{(r/10^i/)-0.25}. After ~ 10* years, the cluster 
settles into its long-term phase of slower internal dynamical 
evolution and stellar evaporation. 

In addition, the observational evide nce so far 
iMackev fc Gilmore |[200^ : lBastian et al. II2008D . though ad- 
mittedly still sketchy, indicates that the extremely young 
clusters show a smaller spread of core radii than the older 
ones. It is therefore tempting to see the shape of the present- 
day SDF (see again Figure lll|l as the result of cluster-to- 
cluster differences in the star formation efficiency, starting 
from a population of protocluster cores that began with 



* A proviso to this statement is that will grow slowly after core 
collapse, which typically occurs after about 20 relaxation times; 
see, e.g., Trenti et al. 2007 among many modelling papers. 
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rather similar sizes. The single-peaked but slightly asym- 
metric shape of the size distribution is strongly reminiscent 
of a normal probability distribution that has passed through 
a nonlinear transformation. In this case, the transformation 
is the conversion of a given SFE to a radial expansion factor. 

In this view, the peak frequency at ~ 2 — 3 pc would 
simply represent protoclusters which experienced the most 
common (i.e. most probable) SFE. Smaller values of the ini- 
tial SFE would lead to clusters with a larger present-day size, 
and the very largest clusters would be ones with SFEs not 
much above the minimum level required to remain bound. 
At the opposite extreme, protoclusters with unusually high 
SFE would suffer only very small expansion and end up at 
the minimum observed size of ~ 1 pc. Said differently, the 
"rise point" of the SDF at ~ f pc should therefore be 
close to the typical initial size of the protoclusters. 

The actual mean SFE within any one protocluster 
should in principle be determined by several factors (e.g. lo- 
cal temperature, pressure, turbulence, degree of initial mass 
segregation), and thus could vary stochastically from one 
core to another. Ho wever, the mean SFE is expected (cf. 
iLada fc Lad"a1l2003l ) to be near 30% for dense star-forming 
regions that will give rise to massive clusters. This mean 
SFE already permits an independent estimate of the amount 
of radial expansion to be expected, from analytical energy 
arguments: in the limit of slow, adiabatic expansion dur- 
ing gas ex pulsion, the product (cluster mass x radius) stays 
constant (|Hills IITqSOI ). which yields an expansion factor of 
{SFE)~^ ~ 3.3. This suggests in turn that the initial sizes 
of the protoclusters should be near ~ 0.8 pc to produce a 
present-day median around 2.5 pc, consistent with the ob- 
servations cited above; 

The models by iBaumgardt fc Kroupa I (|2007l ) (BK07) 
allow us to explore some simple simulations of the SDF 
a bit further along the lines outlined above. BK07 use 
three critical factors determining the final expansion ratio 
ER = rhifinal) /rh{initial): (a) the SFE in the protoclus- 
ter, (b) the initial cluster size relative to its tidal radius 
(rh/rt), and (c) the ratio of the "mass loss time" tm over 
which gas is expelled relative to the internal crossing time 
tc- As argued above, we expect the protocluster core to be 
typically ~ 0.5—1 pc in size. This level is far smaller than the 
typical tidal radius of > 50 pc for a massive globular clus- 
ter, so for our purposes we assume a small (rh/rt) < 0.03 
(see Figure 4 of BK07). In addition, the very large amount 
of dense gas present in the core during star formation will 
be capable of absorbing and thermalizing the SNe ejecta 
and stellar winds, prevent ing instantaneous mass loss (see 
also iBailin fc Harris 1 12009| ) . In these massive protoclusters 
as well, the c rossing time tc ^ 10^ y is less than the gas ex - 
pulsion time (|Goodwin II 19971 : IBaumgardt fc Mieske II2008I ). 
We therefore use the BK07 models for the two representa- 
tive cases Tm = tc (moderately rapid but not instant gas 
expulsion) and tm — lOtc (slow, near-adiabatic gas expul- 
sion) . 

It should be noted that we do not expect these setup 
arguments to carry over identically for low-mass clusters. 
For these, the smaller tidal radius will mean more rapid 
escape of stars during the expansion phase as well as more 
rapid expulsion of gas (|Baumgardt fc Mieske |[2008l ) , both of 
which put the early evolution into a different regime of the 
model parameter space. At masses well below the normal 
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Figure 18. The scale size distributions for the globular clus- 
ters in M104 (solid histogram) and the Milky Way (dashed his- 
togram), repeated from Figure [TT] In the left panel, the solid red 
line shows a simulated distribution for model clusters with slow 
gas expulsion and a range of star formation efRciencies averaging 
SFE = 0.28, cr^; = ±0.11 as described in the text. In the right 
panel the red line shows a simulated distribution for rapid gas 
loss and a mean SFE = 0.30, ctb = 0.07. 



GC range, a far higher fraction of the protoclusters will not 
survive this initial mass-loss stage. 

With these assumptions about the tidal limit and mass 
loss time, the dispersion in the resulting cluster sizes starting 
from a fixed initial core size will be due simply to a disper- 
sion in the SFE's from one cluster to another. To model 
the SFE in a simple way, we assume that it follows a Gaus- 
sian distribution with a mean of Eq and a standard devia- 
tion of (T_B. In general, this approach r esemb les the grid of 
simulations bv IBaumgardt fc Mieske"! (|2008l ) drawing from 
the same models. However, their work was directed towards 
exploring the trends of mean cluster mass and radius as 
functions of galactocentric distance. Here, we concentrate 
on attempting to reproduce the detailed shape of the SDF 
itself. 

Using the BK07 model grid, we find that the relation 
between the formation efficiency and the expansion ratio ER 
can be well approximated by ER ~ 0.55(SFE )~^'^ for th e 
rapid case tm = tc- As already noted above ('cf. lHills1ll980l ). 
ER — {SFE)~^ for the slow-expulsion case tm = lOtc. We 
then run Monte Carlo realizations with these assumptions, 
and vary Eo,oe, and the initial cluster size rh{init) to find 
simulations that match the real distributions. A final bit 
of input is to convert the initial cluster radius rh (three- 
dimensional) to the projected rh that we observe with the 
correction factor rh{proj) = 0.73r;i(3D). 

Figure fTSl shows the results of two sample runs, matched 
to the data shown previously for M104 and the Milky Way. 
Each simulation is the result of generating 10^ clusters at 
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random following the prescriptions described above. The left 
panel shows the "slow expulsion" case for the parametric 
values Eo = 0.28, (Je ~ 0.11, and rh{init) = 1.0 pc. The 
right panel shows the "rapid expulsion" case for Eq = 0.30, 
cr_B = 0.07, and rhiinit) = 0.9 pc. Both the simulations and 
the Milky Way distribution have been normalized to the 
total population of the M104 system. Both models match the 
real data encouragingly well. The three important features 
are the rapid ramp-up in numbers of clusters starting at 
rh ~ 1 pc; the moderately broad peak near 2.5 pc; and the 
long but low-amplitude tail extending to much larger radii. 

A limitation of the present discussion is obviously that 
we have not properly included the long-term effects of dy- 
namical evolution on the SDF. These effects would progres- 
sively trim the SDF that emerged shortly after the early gas- 
loss phase, gradually removing low-mass or very low-density 
objects. For this reason, our SDF models which predict a 
few more clusters at larger radii than in the present-day 
data are probably not a cause for serious worry, because 
these large-radius, low-density clusters are among the ones 
that would preferentially get destroyed by dynamical evo- 
lution over the subsequent Hubble time. The high— tail 
predicted by the models should therefore be only an upper 
limit to the present-day distribution. 

iDa Costa et al. I l|2009f ) have collected the recent obser- 
vations for the SDF particularly in dwarf galaxies, and de- 
velop intriguing evidence that the SDF in total may be bi- 
modal. The "normal" mode peaked at — 2.5 — 3.0 pc is 
the more well populated of the two, but there is a second 
mode peaked near — 8 pc. There are no obvious obser- 
vational selection effects that would bias discoveries (or size 
measurements) against cluster in the intermediate-size range 
around 5 pc between the two modes, so the reality of the bi- 
modality must be taken seriously. Da Costa et al. propose 
that these may belong to a second mode of cluster formation 
more prevalent in the weaker potential wells of the dwarfs. 

Figure [19] shows a match of our model to the old- 
GC sample compiled from several dwarf galaxies by 
iGeorgiev et al. I (|2009l ) described above, where we attempt 
to match the bimodal distribution noted by Da Costa et al. 
The normal mode (solid line) is the fast-expulsion model for 
-Eo = 0.32, ffE = 0.08, and rhiinit) = 1.15 pc, while the 
"extended" upper mode (dashed line) assumes Eq — 0.30, 
as = 0.03, and rhiinit) = 2.4 pc. Reproducing the upper 
mode requires a distinctly larger initial protocluster radius, 
though interestingly, it needs to be larger by only a factor 
of two. The formal value of as is quite a bit lower than for 
the normal mode, but once again it is less clear how sig- 
nificant this is, because clusters that ended up with very 
large radii rn > 15 pc might not have survived except (as 
also discussed by Da Costa et al.) under the most favorable 
circumstances. Minor variations around these combinations 
of parameters can be found that give similar fits, and the 
present discussion should be taken only as illustrative. 

Lastly, Figure [20] shows the sensitivity of the model 
SDF to changes in both the basic SFE and the initial core 
size rhiinit). The same model fit as in Figure [TBb for rela- 
tively rapid expulsion (solid line) is shown (£"0 = 0.30, cje = 
0.07, rh{init) = 0.9) along with the Milky Way cluster data. 
The dashed and dotted lines in the left panel show the 
same model but now with only the mean SFE changed 
(-Eo = 0.35, 0.25). The right panel shows models for changes 
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Figure 19. The scale size distribu tions for the sara ple of globu- 
lar clusters in dwarf galaxies, from iGeoreiev et al. I lioOft ) as de- 
scribed in the text. The solid (red) line shows a simulated model 
distribution for the rapid gas expulsion case with Eq = 0.32, 
CFE = 0.08, rhiinit) = 1.15 pc. The dashed line shows a com- 
parable simulation but now with -Eo = 0.30, cte = 0.03, and 
rhiinit) = 2.4 pc. 



in only the initial radius (rhiinit) = 0.7, 1.1). These four 
outlying models clearly fail to match the data. The implica- 
tions are that the average starting conditions are constrained 
to rather narrow ranges, ±0.05 or less in Eq and ±0.1 pc or 
less in rh{init). To a large extent, variations in rhiinit) can 
be traded off with variations in -Eq to produce the final range 
of observed sizes. 

In rougher terms, the appropriate ranges for the initial 
conditions can be understood physically as follows. If Eo 
falls below 0.20 — 0.25, then very few clusters survive the 
gas expulsion phase at all. At the opposite end, values of Eo 
higher than about 0.4 mean that most clusters experience 
too little expansion to reproduce the peak at 2.5 pc or to fill 
up the high— tail. For the initial size, values of rhiinit) 
bigger than ~ 1 pc or less than ~ 0.7 pc drive the present- 
day peak well above or below the characteristic 2.5-pc level 
that we need, and also fail to match the observed absolute 
range in the SDF. Finally, the third fitting parameter cr_B is 
essentially used to fine-tune the SDF peak and dispersion, 
once Eo and rhiinit) have been put into the right range. 

These sample realizations are undoubtedly oversimpli- 
fied, and do not by any means represent a thorough ex- 
ploration of the parameter space of the models. This ap- 
proach also ignores other potentially interesting effects such 
as external interactions with other clouds; mergers of young 
clusters; dumpiness and substructure within a protoclus- 
ter; or primordial mass segregatio n. All of these could affect 
the early structural evolution (e.g.lScheepmaker et ah II2OO7I : 
iBaumgardt et al. II2OO8I : [Marks et al. Il2008h . 

Our main point, however, is that it seems possible to un- 
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Figure 20. Left panel: The solid histogram shows the size dis- 
tributions for the Milky Way globular clusters. The red solid line 
is the "best-fit" model for rapid gas expulsion (Eq = 0.30, c^; = 
0.07, Thiinit) = 0.9 pc) while the dashed and dotted lines show 
models for Eq = 0.35 and 0.25. Right panel: Same as in previous 
panel, but here the red dashed and dotted lines show models for 
rh{0) = rfi{init) = 0.7 and 1.1 pc. 



derstand the key features of the globular cluster size distri- 
bution with a relatively simple set of assumptions, and with 
quite plausible fiducial values for the range of star formation 
efficiencies and initial protocluster sizes. Various arguments 
in the literature already favo r a mean SFE of ar ound 30 
percent for star clusters (e .g. iLada &: Ladall2003l ) and an 
initial size less than 1 pc (iBastian et al. I 20081 ). However, 
there are far fewer avenues to quantitative estimates for the 
range of SFE's and initial scale sizes that typified the for- 
mation regions of massive star clusters. The detailed shape 
of the present-day size distribution appears to be one such 
method. 



5 BIMODALITY AND THE 

MASS/METALLICITY RELATION 

M104 was one of the first galaxies in which the intrigu- 
ing correlation between mean luminosity and color along 
the blue sequence was found (Paper I). This trend acts in 
the sense that the most luminous clusters become slightly 
but progressively more metal-rich, and thus corresponds to 
a mass/metallicity relation (MMR). The original discov- 
ery papers used sample s of data from large elliptical galax- 
ies (IHarris et al. I l2006l : IStrader et aTl l2006l : iMieske et aTl 
I2OO6I ). but the M104 data indicated that it might extend to 
disk systems as well. The first round of papers gave different 
results for the detailed shape of the MMR, b ut more recent 
data analysis and discussioii s l|Harris |[2009al lbl: iPeng et al. I 
I2OO9I : ICockcroft et al. Ilioogj ) indicate more of a consensus 



Figure 21. Histogram of {B — R) colors for the M104 globular 
clusters in the magnitude interval 20.5 < i? ^ 23.5 (solid points 
with errorbars) . The solid curve shows the best- fit bimodal Gaus- 
sian, with the two individual modes shown as dashed lines. 



emerging around the view that the blue sequence is nearly 
vertical for lower luminosities (L < 5 x 10'' Lq) but gradually 
slants more toward the red going up to higher luminosities. 
Equally intriguingly, no such MMR seems to affect the red 
sequence, which keeps the same mean metallicity at all lu- 
minosities. 

The most effective physical interpretation so far is based 
on some fo rm of self-enrichment during a cluster's forma - 
tion stage l|Strader fc Smith Il200i : iBailin fc Harris 1 120091 '). 
A protocluster of mass > lO^M© or more, within a core of 
~ 1 pc, will be able to hold back enough of its first round 
of SN-enriched material to enrich the still-forming low-mass 
stars, thus giving the entire cluster a higher mean metal- 
licity than the pre-enriched level it started with. This extra 
self-enrichment will be much less noticeable along the red se- 
quence because its protocluster gas is an order of magnitude 
higher in heavy-element abundance than the blue sequence. 

We have used our new photometry of the Ml 04 clus- 
ters to re-investigate the existence of an MMR within this 
massive disk galaxy. As noted above, all the BVR total mag- 
nitudes are individually corrected for the scale sizes of the 
clusters and thus are free of any aperture-size or PSF-fitting 
effects that might depend on luminosity. The first test is 
to measure any correlation of mean color with luminosity 
along both sequences. To do this, we use the {B — R) data 
as the most metallicity-sensitive of the three possible color 
indices we could define, and then divide the data into half- 
magnitude bins by R magnitude. The color distribution is 
then put into t he fitting routine RMIX (|Wehner et al. II2OO8I : 
[Harris ll2009al ) and the best-fit bimodal Gaussian distribu- 
tions are found in each independent interval. In this way, 
we do not assume any particular form for the MMR along 
either sequence. 
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Table 2. Bimodal Fits to the Color Distributions 



R Range 


{R) 


n 




M2 


Blue Fraction 


22.5-23.5 


22.92 


95 


1.118 ± 0.013 


1.479 ±0.015 


0.434 ± 0.054 


22.0-22.5 


22.21 


103 


1.121 ± 0.029 


1.448 ±0.031 


0.463 ±0.102 


21.5-22.0 


21.74 


131 


1.162 ± 0.016 


1.452 ±0.022 


0.489 ± 0.070 


21.0-21.5 


21.29 


116 


1.169 ± 0.022 


1.449 ±0.035 


0.507 ±0.104 


20.5-21.0 


20.77 


82 


1.175 ± 0.013 


1.494 ±0.017 


0.549 ± 0.059 


20.0-20.5 


20.28 


55 


1.214 ±0.017 


1.470 ±0.027 


0.548 ± 0.087 


19.5-20.0 


19.78 


36 


1.178 ± 0.018 


1.499 ±0.028 


0.561 ±0.088 


19.0-19.5 


19.32 


17 


1.258 ± 0.032 


1.518 ±0.047 


0.522 ±0.153 


18.0-19.0 


18.68 


16 


1.278 ± 0.022 


1.520 ±0.061 


0.733 ±0.124 



The mean points derived from these objective fits are 
listed in Tableland plotted in Figure[22l In principle, RMIX 
can solve for five free parameters: the mean colors /ii,/i2 of 
the blue and red modes; their Gaussian dispersions ai,a2; 
and the proportion pi (or p2 — 1 — pi) that the blue (red) 
mode makes up of the total population. In practice, if the 
number of datapoints in the bin is less than about 50, the 
solutions need to be partially constrained for convergence; 
here, we choose in such cases to fix the dispersions (Ti,2 be- 
cause these do not change noticeably along the sequence 
and we are primarily interested in tracing the mean colors 
themselves. By using the total data over 20.5 < R < 22.0 
we find (cri) = 0.070 ± 0.006 and (era) = 0.096 ± 0.010. 
The color distribution in {B — R) for this galaxy defines a 
remarkably clear bimodal histogram (Figure I21|l , and the 
double- Gaussian model provides an accurate fit. 

As can be seen from Table [2] the blue and red clusters 
make up nearly equal proportions of the total population; 
the pi (blue) ratio increases steadily toward higher luminos- 
ity and the blue sequence reaches a bit higher at the top 
end. As expected from both the previous literature and the 
self-enrichment model, the red sequence does not show a 
significant change in mean color with luminosity (that is, it 
has a "zero" MMR). By contrast, the blue sequence shows a 
clear slope toward the red that is nearly linear in form. An 
unweighted linear fit to the mean points in Table [2] gives a 
highly significant slope A{B - R)/AMr = -0.037 ± 0.005, 
almost identical with what was derived in Paper I. In terms 
of heavy-element abundance Z this slope corresponds to a 
scaling Z ~ L°-^^±OM ^^^^ below). 

The slope of the linear fit becomes progressively less 
significant as the higher-luminosity bins are removed. For 
example, a direct fit of {B — R) versus R for the clusters with 
Mr > -8 and {B - R)o < 1.20 gives A(B - R)/AMr = 
—0.013 ±0.014. Thus we cannot rule out the possibility that 
the blue sequence may become more nearly vertical at lower 
luminosities. Nevertheless, we confirm the basic trend found 
in Paper I that the blue sequence does show an MMR. In 
addition, the metallicity scaling agrees extremely well with 
the mean slopes Z ~ L"'^ found for several giant E ga laxies 
in the recent studies bylHarris I (|2009al ) , [Harris 1 (|2009bh ,and 
ICockcroft et al. I l|2009l ). and gives some additional support 
to the idea that the MMR may be a near- "universal" phe- 
nomenon which requires a broad-based physical explanation 
independent of galaxy type. 

As a more direct comparison between theory and data, 
we show in Figure [22b the same mean points from the color- 
magnitude diagram, but now superimposed on a simula- 
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Figure 22. Left panel: Color-magnitude data for the globular 
clusters in M104, from Figure 1101 and transformed to absolute 
magnitude Mr versus intrinsic color {B — R)o- Along each se- 
quence the mean points in half-magnitude intervals (Table [2]| are 
plotted as the connected large dots. The large star at top marks 
the UCD. Right panel: A simulated set of data from a bimodal 
metallicity distribution, as described in the text. At high mass 
the clusters are noticeably affected by self-enrichment. 



tion drawn from the iBaihn fc Harris I (j2009t ) theory. In this 
model, the clusters are assumed to form from dense proto- 
clusters of mass Mo, initial size ro, and with star forma- 
tion efficiency Eq. The protocluster gas has a pre-enriched 
metallicity level [m/H]o, which is then enriched further by 
some fraction of the first SNe that go off in the emerging 
cluster. After its formation stage, the early gas expulsion 
leaves the cluster with a mass Mo x Eq. Then to take ac- 
count of the longer-term dynamical mass loss that follows, 
we use the conventional expression for a roughly constant 
mass loss rate applying to two-body relaxation in tidally lim- 
ited clusters, M(t) ~ Mq — /Jetii (e.g. iBaumgardt fc Makino I 
l2003l : iMcLaughhn fc Fall II2008I ). The actual mass loss rate 
Hev is expected to be a fun c tion o f cluster density; for ex- 
ample, iMcLaughhn fc Fall I l|2008h derive fiev = IWOpl^^ 




Globular Clusters in MIO4 19 



Mq Gyr^^ where ph is in Mq pc~'^. The empirical evi- 
dence shows that the characteristic density ph in turn in- 
creases systematically with mass roughly as pu ~ M, though 
it shows large cluster-to-cluster scatter (see Figure la of 
McLaughlin and Fall, for example). A fully detailed sim- 
ulation of these effects is beyond the scope of our discus- 
sion, but as a first-order representation of the mean trend, 
we use ^ 5OOO(Af/lO^M0)^/2 M© Gyi~^, which rea- 
sonably re produces the mass loss r a tes in the N-body simu- 
lations by Baumgardt fc Makinol ||2003) and more recent 
work of Kruiissen fc Portegies Zwart I 1 20091) and Hurley 
(private communication). To set the other parameters in the 
enrichment model, for all clusters we use Eo — 0.3, ro — 1 
pc (Section 4 above). We also use a supernova conversion 
efficiency fsN = 0.3 (that is, 30% of the SNe in the cluster 
formation sequence happen soon enough in the burst to help 
enrich the remaining gas). 

The initial (pre-enriched) values for the metallicities of 
each sequence are chosen to match the two observed se- 
quences, at {B - R)o = 1.09 (blue) and 1.39 (red). (Note 
here that these are the dereddened colors, which we use for 
the simulations and plot in Figure (22] The data in Table 
2 are the directly observed colors with no reddening re- 
moved.) Transformation between metallicity and color is de- 
term i ned by (B - R)o = 27fFe/H] + 1.52 (|Barmbv et al. I 
l2000l : ICantiello et al. l |2007|). We assume a mass-to-light ra- 
tio {M/L)ii — 2 and convert to absolute magnitude with 
Mr — 4.29 — 2.51og(L/L0). Finally, the sequences are as- 
sumed to have intrinsic dispersions in metallicity of ±0.26 
dex (blue) and 0.36 dex (red). 

With this input, we then construct Monte Carlo realiza- 
tions of the entire GCS, a typical example of which is shown 
in Figure 122b . The initial mass distribution is drawn ran- 
domly from a simple power law form dN/dMo ~ M^'^, and 
the total number of simulated clusters is chosen to match 
the observed numbers in M104 for Mr < —8. Because of 
our overly simplistic initial mass function for the input clus- 
ters, the simulation overestimates the numbers both at very 
high and low luminosity; however, the important feature is 
the matchup in the mean position of each sequence. 

We find that the self-enrichment model can reproduce 
the red and blue sequences basically well. On the blue 
sequence, the "cut-on" point where self-enrichment starts 
clearly to dominate over the pre-enriched level is at Mr ~ 
— 10 (M ~ 10^ Mq), above which the sequence starts to slant 
more strongly to the red. At very high luminosity, both se- 
quences converge towards an enrichment level [m/H] ~ —0.4 
or {B — R)o ^ 1.4, very near where the M104 UCD sits in 
the color-magnitude diagram. 

A second-order concern is the discrepancy between 
model and data for the low- luminosity range Mr > — 8. The 
real clusters in this range appear to be slightly more enriched 
(redder) than their simulated counterparts. The model, as it 
stands, firmly predicts a "zero MMR" at present-day masses 
below about 10^ M©; even in their more massive protocluster 
state, such clusters are simply not massive or dense enough 
to hold back significant amounts of the enriched SNe ejecta. 
In other words, the model provides a good quantitative ex- 
planation for the MMR at the upper end of the sequences, 
but cannot explain any MMR that might continue to lower 
masses. Thus if the MMR genuinely extends to lower lu- 
minosities (which we take our data to hint at, though not 



conclusively), either the model is incomplete or another ex- 
planation must exist. 

The main free parameters in the model are the SFE Eq, 
initial size ro, and SN conversion efficiency fsN- In princi- 
ple, we could adjust these to force the blue cut-on point 
to lower luminosities: for example, making the SFE mass- 
dependent such that Eo is higher at lower masses, or mak- 
ing ro smaller for the lower-mass clusters, would raise the 
resulting enrichment at lower mass. We could also make fsN 
arbitrarily higher only for the lower-mass clusters, but it is 
not clear what might produce such a counterintuitive ef- 
fect. Decreasing ro to 0.5 pc or increasing fsN to 1.0 have 
about the same effect of lowering the cut-on point by less 
than a factor of two, and combining both effects would re- 
duce it to 4 X 10^ Mq. However, changes of this type away 
from their nominal values of Eo — 0.3 and ro — 1 pc would 
dramatically alter the resulting SDF and are thus ruled out 
by the analysis presented in Section 4 above. In addition, 
they would increase the self-enrichment at higher luminos- 
ity even further, and destroy the matchup with the real data 
for M > 10^ Mq. 

For the present time, we see no reason to discard the 
basic approach of the self-enrichment model, which is quite 
effective at matching the central features of the MMR par- 
ticularly for the high-luminosity clusters. The present study 
should be viewed as only an indication that further inter- 
esting tests may rely on the detailed features of the GC 
metallicity distribution over the entire run of luminosities. 



6 SUMMARY 

In this study we have used a mosaic of EST/ ACS images of 
the Sombrero galaxy (M104) to make new measurements of 
its large globular cluster population. We obtain new BVR 
photometry and effective-radius measurements for a total of 
652 clusters. Our findings include the following: 

(1) The individual cluster profiles have been fitted with King 
and Wilson dynamical models to derive their effective (half- 
light) radii r^. A series of tests of the data with different 
fitting models and different definitions of the stellar point- 
spread function show that the internal precision of the half- 
light radii is typically ±0.4 pc. 

(2) The distribution function of the effective radii (SDF) 
has a peak near r^ — 2.4 pc and an asymmetric tail to 
larger radii, which to first order is remarkably similar to the 
SDFs found for other galaxies, including the Milky Way, 
giant E galaxies, and several kinds of dwarf galaxies. We 
find second-order trends with cluster metallicity and spatial 
location (the clusters are slighter bigger on average at lower 
metallicity, or larger galactocentric distance). Both of these 
trends also resemble what has been found in other galaxies 
in several other recent studies. However, closer comparisons 
with other data in the literature show that the metallicity- 
based size difference may not be as large for GCs in disk 
galaxies as in E galaxies. 

(3) The M104 GCs define very much the same regions in 
the structural Fundamental Plane as in other, more well re- 
solved systems including the standard Milky Way. We use 
the half-light surface brightness Iv,h and binding energy Ei, 
to demonstrate that clusters become preferentially more ex- 
tended at larger galactocentric distances or low luminosities. 
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As long as the cluster half-light radii can be measured ac- 
curately, these quantities can be used in still more distant 
systems such as the giant ellipticals in Virgo and elsewhere. 

(4) We have explored a simple framework for the phys- 
ical origin of the globular cluster scale-size distribution. 
We assume that GCs start as dense, massive protoclusters 
which form stars at a certain efficiency and then expand 
to radii near their present-day sizes during a rapid initial 
st age of mass loss and residual g as expulsion. Using models 
bv lBaumgardt fc Kroupa I (|2007l ). we find that the observed 
SDF can be successfully and closely matched if we assume 
that (a) the protoclusters began with scale sizes ~ 0.8 — 0.9 
pc; (b) the mean star formation efficiency was ~ 0.30 but 
with stochastic cluster-to-cluster variations on the order of 
±0.08 — 0.15; and (c) the gas expulsion time was at least as 
long as the internal crossing time. 

(5) The color-magnitude distribution for the GCs shows a 
clearly defined, classic bimodal form with nearly equal num- 
bers of metal-poor and metal-rich clusters. Detailed bimodal 
fitting of the {B — R) colors shows that the blue, metal- 
poor sequence exhibits a well determined mass/metallicity 
relation (MMR), becoming slightly but steadily redder to- 
wards higher luminosity. The red sequence is vertical over 
the entire luminosity ra nge. Detailed compariso ns with the 
self-enrichment model of iBailin fc Harrlsl (|2009l ) show that 
a close first-order match to the data can be obtained with 
the same choice of star formation efficiency and protocluster 
size that are required to successfully model the SDF. How- 
ever, a potential problem for the model may lie in whether 
or not the MMR extends further down the blue sequence 
below the point where self-enrichment can be expected to 
work. This will be an intriguing area for future work. 
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